Strain, Magnetic Anisotropy, and Anisotropic Magnetoresistance in (Ga,Mn)As on 
High-Index Substrates: Application to (113) A-Oriented Layers 
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Based on a detailed theoretical examination of the lattice distortion in high-index epilayers in 
terms of continuum mechanics, expressions are deduced that allow the calculation and experimental 
determination of the strain tensor for (/i/i/)-oriented (Ga,Mn)As layers. Analytical expressions are 
derived for the strain-dependent free-energy density and for the resistivity tensor for monoclinic 
and orthorhombic crystal symmetry, phenomenologically describing the magnetic anisotropy and 
anisotropic magnetoresistance by appropriate anisotropy and resistivity parameters, respectively. 
Applying the results to (113)A orientation with monoclinic crystal symmetry, the expressions are 
used to determine the strain tensor and the shear angle of a series of (113)A-oriented (Ga,Mn)As 
layers by high-resolution x-ray diffraction and to probe the magnetic anisotropy and anisotropic 
magnetoresistance at 4.2 K by means of angle-dependent magnetotransport. Whereas the transverse 
resistivity parameters are nearly unaffected by the magnetic field, the parameters describing the 
longitudinal resistivity are strongly field dependent. 
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I. INTRODUCTION 

For the past two decades, dilute magnetic semicon- 
ductors, and in particular ferromagnetic semiconductors, 
have been attracting considerable attention due to their 
exceptional physical properties as well as their potential 
applicability in information technology. Ferromagnetism 
mediated by delocalized p-type carriers- could be imple- 
mented in the standard semiconductor GaAs by incorpo- 
rating magnetic Mn atoms on Ga sites, cf. Ref. 2; and 
references therein. Since the highest Curie temperature 
reported so far is 185 Kf^ the application of (Ga,Mn)As in 
electronic devices operating at room temperature seems 
to be doubtful. Nevertheless, due to its specific electronic 
and magnetic properies, (Ga,Mn)As represents an ideal 
test system for future spintronic applications. 

Magnetic anisotropy (MA) and anisotropic magne- 
toresistance (AMR) are well established key features of 
(Ga,Mn)As, largely arising from spin-orbit coupling in 
the valence band. Consequently, both, MA and AMR 
strongly depend on crystal symmetry and strain. On the 
one hand, this dependence defines the magnetic hard and 
easy sx.esr'^ controlling the orientation of the magne- 
tization and thus the electrical resistivity. On the other 
hand, it offers a method to intentionally manipulate the 
magnetic properties, e.g. by applying external strain via 
piezoelectric actuators.^ In any case, a quantitative de- 
scription of the relation between crystal structure and 
MA and AMR, ideally by means of analytical expres- 
sions, is imperative. 

So far, most of the published work focuses on 



(Ga,Mn)As grown on (OOl)-oriented substrates while 
only few publications report on (Ga,Mn)As grown on 
high-index substrates such as (113) and (114)^— where 
the description of the MA and the AMR is more compli- 
cated due to the reduced crystal symmetry. Potential ap- 
plications of high-index (Ga,Mn)As are ridge structures 
with (113)A sidewalls and (001) top layers, which could 
be used, e.g. for memory devices exploiting different co- 
ercitive fields of the sidewalls and top layers. Yet, a co- 
herent theoretical description of the MA and AMR for 
layers on high-index substrates, which takes into account 
the symmetry of the strain tensor e, is still missing. 

In this work, we present a concise phenomenological 
description of the MA and AMR for (ft,ft,^)-oriented fer- 
romagnetic layers. To this end, we describe the struc- 
tural properties of epitaxial high-index layers in terms of 
continuum mechanics and apply the theoretical results 
to the case of (113)A orientation. High-resolution x-ray 
diffraction (HRXRD) measurements were performed on 
a series of (113)A-oriented (Ga,Mn)As layers to quan- 
titatively determine the epitaxial strain and symme- 
try of the layers. The MA and AMR were investi- 
gated by angle-dependent magnetotransport measure- 
ments. In agreement with previous studies on (113)A- 
oriented (Ga,Mn)As)^""i2 we observe a uniaxial MA along 
the [001] direction, which can now be explained in the 
light of our theoretical model. Whereas the transverse 
resistivity parameters turn out to be nearly constant, a 
systematic dependence of the longitudinal resistivity pa- 
rameters on the strength of the external magnetic field is 
found. 
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The paper is organized as follows: In Sec. Ill Al we 
calculate the strain tensor E for (/i/iZ)-oriented layers 
and provide formulae that allow experimental access to 
e. The results are applied to (113)-oriented (Ga,Mn)As 
samples where the Bravais lattice is base-centered mon- 
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oclinic. An extension of the theory to partially relaxed 
layers is given in Appendix |^ In Sec. IIIBi we present 
a phenomenological expression for the MA taking into 
account the specific form of the strain tensor for [hhl)- 
oriented layers. In Sec. Ill C[ we derive expressions for the 
longitudinal and transverse resistivities piong and ptrans, 
respectively, which apply to monoclinic crystal symme- 
try and current direction along [332] . The resistivity ten- 
sors for monoclinic and orthorhombic symmetry can be 
found in Appendix [Bj They are required in the deriva- 
tion of piong and ptians for arbitrary current directions. In 
Sec. nil AI and Sec. IIIIBI the experimental results of the 
HRXRD and magnetotransport studies are presented. 

II. THEORETICAL CONSIDERATIONS 

Crystal symmetry and epitaxial strain strongly influ- 
ence the MA and the AMR. Therefore, we start with a 
detailed theoretical examination of the lattice distortion 
of {hhl) layers in terms of continuum mechanics and ap- 
ply the general results to (113)-oriented layers (Sec. Ill At . 
Based on these results, we present phenomenological ex- 
pressions for the MA in strained high-index ferromag- 
netic layers (Sec. IIIBI) . We present the resistivity ten- 
sors for monoclinic and orthorhombic symmetry (Ap- 
pendix [BJ and calculate piong and ptrans for monoclinic 
symmetry and current direction along [332] fSec. Ill C|) . 

A. Structural Properties 



U strained layer 




FIG. 1: (color online) The points S, O, and R are located 
within the interface between substrate and layer, n, p, m, 
and m' are unit vectors, n and p are normal and parallel 
to the interface, respectively; m and m' axe perpendicular 
to the lattice plane m — (hkl) and m' = {h'k'l'). They lie 
within the plane spanned by n and p and are symmetrically 
oriented with respect to n (r^ = t^')- 

Before any strain is being applied, the point T has the distance 
X from the interface. The displacement to its final position U 
can be described as a superposition of a hydrostatic compres- 
sion of the (unstrained) layer (T — >■ P) and a displacement 
by xa (P — >■ U), cf. Eq. U]). The hydrostatically compressed 
layer has the same lattice plane spacing as the substrate and 
therefore Eq. ((S| and ^ can be derived considering only the 
displacement described by a. 



1. Distortion of Arbitrarily Oriented Epitaxial Layers 

The strain in an epitaxial layer can unambiguously be 
described by the distortion tensor A with components 
Aij ~ du,;/dxj, where u is the mechanical displacement 
field and Xi denote Cartesian coordinates along the cubic 
axes. In order to calculate the distortion of epitaxial 
layers grown on arbitrarily oriented cubic substrates, we 
decompose A according to Hornstra and Bartelai^ 



A = leii -I- a (g) n, 
where ® denotes the dyadic product and 



Eh 



ai 



(1) 



(2) 



is the isotropic strain that compresses (expands) the cu- 
bic unit cell of the layer to the size of the substrate's unit 
cell; as and ai denote the lattice parameters of substrate 
and relaxed layer, respectively, n is the unit vector per- 
pendicular to the interface and a is a vector that repre- 
sents the anisotropic distortion of the layer. Figure [T] 
illustrates the superposition of the displacements de- 
scribed by Eq. ([T]) : If we consider any mathematical point 
in the unstrained layer at the distance x from the inter- 
face, then this point has the distance (1 + eh)x from the 



interface after applying Ish- The displacement from this 
position into its final position is given by ax. 
If the layer is in static equilibrium, there is no stress 
perpendicular to the surface. Applying Hooke's law, this 
constraint leads to the set of linear equations (z = x, y, z) 



CijkiSkinj = ^ (^£hSki + -^{akni + a/rtfc)^ Cijkinj = 0, 

(3) 



jki 



where 



£kl 



Am + Aiu 



(4) 



are the components of the (symmetrized) strain tensor 
and Cijki denote the elastic stiffness constants of the crys- 
tal. With given Cijki and n, Eq. ([3]) can be solved for 
the components a^, which are proportional to Eh- Con- 
sidering Eq. ([T]), the distortion tensor A is then known 
as a linear function of the only remaining parameter eh, 
which can be determined experimentally as follows. 
We label lattice planes by a single index m which stands 
for a Miller index triplet, e.g. m = (hkl), and we some- 
times refer to planes by their normal m. A lattice plane 
m within a distorted layer in general encloses an angle 
At™ with the corresponding lattice plane in the sub- 
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FIG. 2: (color online) (a) Magnitude 6 = |a/eh| and (b) angle 
/3, calculated with various stiffness constantsi^, are plotted as 
a function of the angle a between n and [001]. Only for the 
highly symmetric crystal orientations [001], [111], and [110] a 
aligns with n. 



strate, cf. Fig. [T] Furthermore, the distortion causes 
a relative difference {lS.d/d)jn between the lattice plane 
spacings in layer and substrate. Both, {Ad/d)m and 
Arm, are accessible via HRXRD and are related to £h 



by the expressions 

Ad 
~d 



m ■ a(eh) 

cos T„ 

1 + eh 
m ■ a(eh) . 

— Sm Tm 



1 



(5) 
(6) 



In the remainder of the paper, we focus on the practically 
relevant case of (/i/i/)-oriented substrates. Solving Eq. ([3]) 
for substrate orientations n between [001] and [110], the 
vector a reads as 



sin(a + /3)/^/2 \ 
sin(a + /3)/\/2 eh , 
cos(q! + 13) ) 



(7) 



where 6 denotes the magnitude of the vector —a/sh and 
/3 the angle between n and —a/eh- In Fig. O b and /3 
are plotted for several cubic semiconductors as a func- 
tion of the angle a between n and [001]. The vector a 
always lies within the (110) plane; thus, this plain is a 
symmetry element (mirror plane) for all layers grown on 
(hhl) substrates. Employing Eqs. (dJ, (g]), and ([7]), the 



components of e can be inferred from Fig. [2] using the 
equations 



-yy 



1 — - sin(a + /3) sin a 



Eh [1 — & cos(a + /3) cos a] 
^xx eh 
bsh 



sin(Q! + 13) cos a 



2V2 

+ cos(a + /3) sina 



(8) 

(9) 
(10) 



(11) 



Equations (l8))-(ITTt are essential for the understanding of 
the MA and will be used in the derivation of the free- 
energy density in Sec. Ill Bl 

Figure [5] shows that for crystal facets other than (001), 
(110), and (111), a is not aligned with the surface normal 
n and the layer is therefore sheared towards a direction 
given by the projection of a onto the surface (cf. Fig.[T]). 
Employing Eq. ^ , the shear angle <^ of a n-oriented layer 
towards any direction p _L n, i.e. Tp = 90°, is obtained 
from 



P ■ a(eh) 



1 



(12) 



where eh has to be determined by the procedure de- 
scribed above. can also be measured directly without 
making use of the stiffness constants and the explicit form 
of a (eh), if we choose two lattice planes m and to' with 
p = (m — m')/2sinrm; insertion into Eq. (fT2|) yields 



a{m — m') 
2sinT„i(l + eh) 



Ar„ 



Ar„ 



Eq. m 



2 sm T„ 



(13) 



in agreement with Ref.|Tj. The angles Ar^ and Ar^/ can 
be derived from rocking curves as described in Sec. IIII Al 
Equation is valid if the two lattice planes m and to' 
are equally inclined towards the surface, i.e. if — r,„/, 
cf. Fig.n 

So far, we have restricted our considerations to pseu- 
domorphically grown layers. With minor modifications, 
Eqs. (O and (jH]) can also be applied to partially relaxed 
layers. In Appendix \^ we discuss how the relaxed lat- 
tice constant and the degree of relaxation can be inferred 
from reciprocal space maps (RSMs) for arbitrarily ori- 
ented substrates applying the formalism described above. 



2. Application to (113)-0riented (Ga,Mn)As Layers 

We now apply the general equations derived in the 
preceeding section to the case of (113)-oriented layers. 
For the following calculations, we use the elastic stiffness 
constants of GaAs given in Ref. Tsl neglecting the Mn al- 
loying. This assumption will be justified for (Ga,Mn)As 
layers with Mn concentrations below 5% by the experi- 
mental results presented in Sec. IIII Al From Fig. [2] we 
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obtain for a = 25.2° the values h = 1.81 and /3 = -17.5°. 
Equation ([T]) then yields 

/ 0.095 \ 

a = -1.81 0.095 eh- (14) 
\ 0.991 ) 

With Eqs. (|8|)- (fTT1) we find for the strain tensor 

/ +0.95 -0.05 -0.35 \ 
e = eh -0.05 +0.95 -0.35 . (15) 
\ -0.35 -0.35 -0.62 / 

Figure [3] schematically illustrates the crystal structure of 
the distorted (113) layer. The position of each atom in 
the strained layer was constructed by applying Eq. ([1]) 
together with Eq. (fT4l) : We start by putting a hydrostat- 
ically compressed layer on the substrate. The unit cell of 
the layer in this hypothetical state is identical to the unit 
cell of the substrate, cf. the dashed and solid cubic unit 
cells in Fig. [3] (a). The final position of each atom can 
be found by a displacement along the direction given by 
a, where the magnitude of the displacement is propor- 
tional to the distance x of the atom from the interface 
(cf. Fig.HD. 

The distortion breaks the symmetry of the layer and the 
only remaining symmetry element is a (110) mirror plane; 
hence, the distorted crystal is assigned to the point group 
m (Cs\ This assignment is important for the derivation 
of the resistivity tensor presented in Appendix [BJ For 
clearness, the crystallographic unit cell of the distorted 
layer is depicted in Fig. |3] (b); it becomes evident that 
the corresponding Bravais lattice is base-centered mono- 
clinic. 

In order to connect Eh with the experimentally acces- 
sible quantity (A(i/(i)„i, we apply Eq. ([5]) to the case of 
the symmetric m=(113) reflection (rua = 0°), i.e. we 
consider lattice planes parallel to the surface. Equation 
([5]) simplifies to 

— = -1.723 — — = 1.723-! ^. (16 

a / 113 1 + eii as 

With Eh as obtained from Eq. ([T6| . the shear angle of 
the (113) layer towards [332] is inferred from Eq. ([T2|) . 
reading now as 

= A7332 = -0.544 = 0.544 ~ . (17) 

1 + Eh Os 

Figure [3] (a) tells us that this shear angle (/) of the layer 
as a whole is the same as the angle Ar332 between the 
(332) lattice planes of substrate and layer, respectively. 

For a direct measurement of (/), we use Eq. 
Referring to Fig. |4] and Fig. [U we choose m' — (333) 
and find the corresponding plane m with equal inclina- 
tion with respect to the (113) plane by interpolating be- 
tween the (115) plane and the (004) plane as described 




FIG. 3: (color online) (Ga,Mn)As(113)A layers exhibit a mon- 
oclinic crystal symmetry, (a) The layer as a whole is sheared 
towards [332] by an angle = AT332, leaving the (110) mirror 
plane as the only remaining symmetry element, 
(b) 3-dimensional view of the monoclinic base-centered unit 
cell with respect to the distorted cubic unit cell. For the sake 
of clarity, the second basis atom of the zinc-blende lattice is 
omitted in this sketch. 




FIG. 4: (color online) Normal vectors of the lattice planes con- 
sidered in Eq. (|18p . Referring to Fig. [T] we chose m' = [333] 
and p = [332]. We find a virtual lattice plane m with equal 
inclination with respect to the (113) plane by interpolation 
between the (004) and (115) plane. For the corresponding 
angles we find = T333 = too4 + 0.27(too4 — Tlis)- 
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in Ref. [ij. We obtain for the shear angle 

Aroo4 + 0.27(Atoo4 - ^nu) - AT333 
2sin (r333) 



The parameters bi denote magnetoelastic coupling con- 
stants. Using the trivial identity + + — 1, 
(18) the contribution Bxymxrriy in Eq. (f2T|) can be expressed 



B. Magnetic Anisotropy 

MA is the dependence of a system's free-energy den- 
sity F on the orientation of the magnetization direction 
m — M/M. In the following, we assume the sample to 
consist of a single ferromagnetic domain with a uniform 
magnetization whose magnitude M is assumed to be con- 
stant; we therefore analyze the quantity Fm — F/M. 
For a phenomenological description of the MA in (hhl)- 
orientcd (Ga,Mn)As layers, we expand Fm in powers of 
the components m^, my, and rriz of m along the cu- 
bic axes [100], [010], and [001], respectively. Considering 
terms up to the fourth order in m, the only intrinsic con- 
tribution to Ff^j for an undistorted cubic layer is a cubic 
anisotropy 



F'cub. 
M 



(19) 



due to the crystal symmetry. Extrinsic contributions to 
Fm are the shape anisotropy, caused by the demagnetiza- 
tion field perpendicular to the layer, and a controversially 
discussed^^'^^ uniaxial in-plane contribution 



F^f =Sd(m-n)2 + S„(m-t)2 



(20) 



where n and t denote unit vectors along the surface nor- 
mal and [110], respectively, and is related to the mag- 
netization M by Bd = ij,qM/2. 

For distorted layers, further intrinsic contributions pro- 
portional to the strain components Sij occur. These are 
referred to as magnetoelastic contributions and are pre- 
sented for arbitrarily strained cubic crystals in Ref. [itI . 
Considering intrinsic and extrinsic contributions, the 
free-energy density for (/i/iZ)-oriented layers can be writ- 
ten as 

FM{'m) = const. + B^irri^ + Bxyirixmy 

+ Bxzinix + rny)mz + B;,iml + Bxiim* + my) 
+ Bxyz2mxmyml + B^2j^^(toj. -|- my)mxmymz 
+ BA{m ■ nf + B^{m ■ tf . (21) 

The anisotropy parameters are related to the strain com- 



ponents by 








B,2 


= bi{£xx 


- £zz) 


(22) 


Bxy 


= 2b2£xy 




(23) 


Bxz 


= 2b2£xz 




(24) 


Bx^ 


— b4Exx~ 




(25) 




= bi£zz — 


363 + 264 , , , 
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(26) 


B 2 


— 2b^£xy 


(27) 


Bx2z2 


= ^b^Exz- 




(28) 



in terms of B^{m ■ t) and B^ 



Therefore, we will 



understand it to be contained in the latter terms in the 
analysis of the experimental data in Sec. IIII BI 

In the case of (OOl)-oriented layers, the off-diagonal 
elements of e vanish and Eq. pT|) simplifies to the well 
known expression 

FM{fn) — const, -f {B^2 + Bd)ml + B^4.m\ 

+ Bx^mt + m%) -f \^B^{mx - myf . (29) 

Equation (j2ip in particular explains in a natural way 
the occurrence of the uniaxial anisotropy Bz^mj, along 
[001], which has been introduced ad- hoc in previous 
publications in order to explain the results of angle- 
dependent ferromagnetic resonance®'^ and magnetotrans- 
port studies^iiiS on (113)A-oriented (Ga,Mn)As layers. 
Furthermore, Eqs. (f^ . and account for the 

strain dependence of the parameters -6^2, -6^4, and B^i^^ 
found in a systematic study of (Ga,Mn)As layers grown 
on relaxed (OOl)-oriented (In,Ga)As buffers^: For (001) 
orientation, the relation Sxx ~ —^■^£zz holds and there- 
fore Eqs. (ESI) and ^ read as 5^4 = {0.6b3-0.7bi)ezz~ 
Bcub/2 and B^4 = (O.663 -I- l.Ab4)£zz - Bcuh/2, respec- 
tively. Figure 10 in Ref. |j shows an increase of Bx4 and a 
decrease of B^t with increasing Szz- These findings agree 
with Eqs. ^ and (ESI if 63 > 764/6 and 63, 64 < 0. 



C. Anisotropic Magnetoresistance 

It is well established that (similar to the MA) the 
AMR, described by the resistivity tensor /?(m), is 
strongly affected by the crystal symmetry. In order to 
obtain an analytical expression for p{m), we performed 
a symmetry-based series expansion of the tensor compo- 
nents up to the fourth order in m. For cubic and tetrag- 
onal symmetry, the explicit form of the resistivity tensor 
and a detailed description of its derivation can be found 
in Ref. 18. In this work, we generalize the expressions 
for p{m) to monoclinic and orthorhombic symmetry. As 
shown in Fig. [lU] of Appendix [B] orthorhombic symme- 
try applies to (llO)-oriented substrates. For other crystal 
facets {hhl) with 7^ /i 7^ the crystal exhibits mono- 
clinic symmetry. The explicit forms of the corresponding 
tensors pirn) are given in Appendix |B] 

The AMR is usually probed by measuring the longitu- 
dinal and transverse resistivities piong and ptrans, respec- 
tively, which are related to p{m) by piong = 3^ ■ P ■ 3 
and ptrans = ' P ' 3 ■ The unit vectors j and t point 
along the current direction and the transverse direction, 
respectively. The resistivities also allow experimental ac- 
cess to the MA as shown in Refs. illll In Sec. HlBl 
the magnetic anisotropy parameters introduced in the 
preceding section are derived experimentally by measur- 
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ing the angular dependence of piong and ptrans at various 
fixed magnetic field strengths. 

Now we turn to (113)-oriented layers. According to 
our previous work/* we are referring to the right-handed 
coordinate system (ji, t,n), where j \\ [332], t \\ [110], 
and n \\ [113]; mj, mt, and r7i„ denote projections of m 
along these directions. We calculate piong and ptrans by 
projecting the resistivity tensor in Eqs. (|B4[) along j and 
t, respectively. We find 

Piong = Po + pim'^ + P2m^n + P3"i| + P4m^ + psm^mf^ 
+ pQiuijiTin + pa2m^jmn + posmjiTil (30) 

and 

Ptrans = P6"i„ + PTm^nit + p^m\ 

+ pgrnjUitml + piom^mt + piimtrrin 

+ pi2mtml + pisTUj + pum^j + pi5m|m„ 

+ piemj-ml + pnm'jmtmn, (31) 

where pi are linearly independent resistivity parameters 
related to the expansion coefficients of p. In the limit of 
unstrained layers (cubic symmetry) the expressions for 
Piong and ptrans havc the same form as in Eqs. ((30|) and 
pi|) . however several resistivity parameters become lin- 
early dependent. This is due to the fact that the current 
direction along [332] already breaks the cubic symmetry. 

AMR studies frequently focus on the special case where 
the magnetization lies within the layer plane (m„ — 0). 
Eqs. dSO]) and ^ then simplify to 

Piong = pa + pim^ + psm'^ (32) 

and 

Ptrans = prmjrfit + piom^mt + pismj + pi4m^.(33) 

III. EXPERIMENT 

We apply the theoretical expressions obtained in the 
preceding section to a series of (113)A(Ga,Mn)As lay- 
ers with different Mn concentrations. In Sec. IIII Al the 
hydrostatic strain e-t and thus e as well as the shear an- 
gle of the layers are derived quantitatively from HRXRD 
measurements. In Sec. IIII Al we present angle-dependent 
magnetotransport measurements which are theoretically 
described by the expressions for piong and ptrans given by 
Eqs. (150)) and (|3ip and by the free-energy density given 
by Eq. 

A. HRXRD Measurements 

(Ga,Mn)As layers with manganese contents up to 
5% were grown on (113)A-oriented GaAs substrates by 
low-temperature molecular-beam epitaxy as described 
in Refs. [l^ and [1^ The structural properties of the 



(Ga,Mn)As layers were experimentally investigated by 
HRXRD. We used a Bruker Siemens D5000HR x-ray 
diffractometer operating with the Cu-Kq^ radiation (A = 
0.154 nm). In order to measure the strain eh and the 
thickness of the (Ga,Mn)As(113)A layers, we performed 
UJ-2Q scans from the symmetric (113) reflection. From 
the angular spacing of the layer thickness fringes^^, we 
inferred a layer thickness of 150 nm. Employing Bragg's 
law, the UJ-2Q scan yields {Ad/d)ii3, and via Eq. we 
obtained the hydrostatic strain e^. With that value ofsh, 
we calculated the shear angle (f> of the layer towards [332] 
using Eq. pT)) . In Fig. [5l eh and (j) are plotted against 
the Mn content. 

In order to verify the consistence of the formalism pre- 
sented in Sec. Ill A 11 we measured the shear angle (j) 
of several samples directly by applying Eq. ^TE\\ . We 
inferred the angles AT333, Aroo4, and Arjis from w- 
scans (rocking curves) with opened detector slits at high 
{+) and low (— ) incidence. For an asymmetric reflec- 
tion, as e.g. the (333) reflection, the corresponding 
lattice plane encloses an angle raaa with the surface, 
and the reflex can be measured at two different angles 
^333 ~ ©333 i ''333 with respect to the surface, where 
6333 is the Bragg angle of the (333) reflection—. If the 
peak separation Aa;333 of layer and substrate is measured 
at high and low incidence, Ar333 can be determined via 
Ar333 — (Awjigg — Aa;333)/2. In Fig.O the results for the 
shear angle (j) obtained in this manner are shown in com- 
parison to those derived from the ^-20 scans. The excel- 
lent agreement of the values confirms the consistency of 
the theoretical formalism. Furthermore, it demonstrates 
that the elastic stiffness constants of GaAs are a good 
approximation for those of (Ga,Mn)As within the inves- 
tigated range of Mn concentrations. 




Manganese content x {%) 



FIG. 5: (color online) Shear angle </!> towards the [332] direc- 
tion and hydrostatic strain eh of (113)A-oriented (Ga,Mn)As 
layers plotted as a function of the manganese concentration. 
The circles represent values for Sh and (j) derived from tj-20 
scans using Eqs. (|16|) and (|17|) . respectively. The stars denote 
values for (j> derived from rocking curves applying Eq. (|18p . 
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B. Magnetotransport Measurements 

Most of the samples under study were found to be in- 
sulating at T = 4.2 K and could therefore not be investi- 
gated by magnetotransport. The hole densities and Curie 
temperatures of the three conducting samples were deter- 
mined from high-field magnetotransport measurements 
as described in Ref. [23I . The results are summarized in 
Table H 

TABLE I: Structural and electronic properties of metallic 
samples, studied via angle-dependent magnetotransport, cf. 

secnm 



3:^Mn(%) 


eh(%) 




rc(K) 


4.9 


-0.26 


2.0 


44 


4.2 


-0.23 


2.2 


47 


3.1 


-0.18 


1.5 


38 



In order to investigate the MA and AMR in the 
(113)A-oriented (Ga,Mn)As samples, we performed 
angle-dependent magnetotransport measurements i^i^iii^ 
For this purpose, the samples were patterned into 
0.3 mm-wide Hall-bar structures oriented along [332] 
with Ohmic Au-Pt-Ti contacts and the longitudinal volt- 
age probes separated by 1 mm. The dc-current density 
was 220 Acm~^. The samples were mounted on a rotat- 
able sample holder in a liquid-He-bath cryostat, which 
was placed between the poles of a LakeShore electromag- 
net. With this setup, the magnetic field could be rotated 
arbitrarily with respect to the crystallographic axes of 
the (Ga,Mn)As layer. 

We rotated the external field H at various fixed field 
strengths within the three different crystallographic 
planes depicted in Fig. Inland measured piong and ptrans- 
In the presence of an external magnetic field H, the free- 




FIG. 6: (color online) Configuration (I), (II), and (III) for 
angle-dependent magnetotransport measurements. The ex- 
ternal magnetic field H was rotated within the planes normal 
to n II [113], j II [332], and t \\ [IlO]. 

enthalpy density Gm — G/M instead of the free-energy 
density Fm determines the magnetization orientation. 
We thus write 



Gm{H) = FM{m) - iiom ■ H, 



(34) 



with F]\{{m) from Eq. ([21]) . The orientation m of the 
magnetization at a given external field H can be found 



by minimizing Eq. (|34)) with respect to m. At the 
maximum applied field of iiqH = 0.62 T, the Zeeman 
term fioHm dominates the free-enthalpy density and 
m essentially aligns along H. With decreasing field 
strength however, the MA described by the anisotropy 
parameters in Eq. (|2ip more and more governs the 
motion of the magnetization as H is rotated with 
respect to the sample. 

By fitting Ec^s. ((30)) and (j3ip to our experimental data 
recorded at fioH = 0.62 T, we obtained values for the 
resistivity parameters pi. Using these parameters, we 
simulated the measured angular dependencies of piong 
and ptrans at weaker fields by varying the anisotropy 
parameters until the simulated curves fit the experi- 
ment. Figure [7] exemplarily shows the experimental and 
simulated angular dependencies of piong and ptians for 
the sample with 4.9% Mn. 

The experimentally obtained anisotropy parameters 
are listed in Table HH In agreement with other workj^iiS 
the MA in our samples can be described by the parame- 
ters B^i, B^i, Bz2, Bd, and B^- As expected, the shape 
anisotropy parameter B^ = increases with in- 

creasing Mn concentration. Since both, the Mn concen- 
tration and the epitaxial strain, infiuence the parameters 
B^i , B^4 , and B^2 , it is not possible to infer an unam- 
biguous strain-dependence of these parameters from our 
experiment. Nevertheless, some qualitative conclusions 
can be drawn: 

Because the relations B^-i — B^i ~ 64 [exx — £zz) ~ 
1.57 64 £h > and Eh < hold for all samples, we find a 
negative magnetoelastic coupling parameter 64 in agree- 
ment with the discussion in Sec. Ill Bl Since the param- 
eters Bxz, Bx^z^^ and Bxyz^ are negligible (the influence 
of Bxy is contained in B^^ and B^ due to the linear de- 
pendence mentioned earlier), we are led to the conclusion 
that the coupling parameters 62 and 65 are smaller than 
bl (at least for the samples with xun > 4%). For more 
strongly strained samples or for samples with larger cou- 
pling parameters bi, however, all anisotropy parameters 
given in Eqs. (|22p - ([27)) may play a role. 
In contrast to our previous experimentsj^ii^ the drastic 



TABLE II: Anisotropy parameters of the samples under study 
as obtained from angle-dependent magnetotransport mea- 
surements. 



2;Mn(%) 


B,4(mT) 


B,4(mT) 


B,2(mT) 


Bd(mT) 


Bu(mT) 


4.9 


-46 


-11 


20 


18 


-8 


4.2 


-56 


-25 


20 


12 


-8 


3.1 


-65 


-40 


5 


5 


-5 



change of the longitudinal-resistivity curves, in particular 
those obtained in configuration (I) , upon variation of the 
external field strength (cf. Fig. [5] and [7]) cannot solely be 
explained by MA. A satisfactory agreement between the- 
ory and experiment can only be obtained by allowing for 



[332] [110] [332] [110] [332] [113] [lio; 



u 

Eh 

03 



Eh 

03 




200 300 400 100 200 300 400 -100 -50 50 
angle of rotation (degree) 

FIG. 7: Angular dependence of the longitudinal and transverse resistivity at various external magnetic field strengths rotated 
within three different planes, referred to as configuration (I), (II), and (III), respectively (cf. Fig. The sample contains 
nominally 4.9% Mn and the hydrostatical strain is Sh = —0.26%. The longitudinal-resistivity curves can be simulated using H- 
dependent resistivity parameters, cf. Fig. [S] In particular, for configuration (I) the longitudinal-resistivity parameters strongly 
vary with the field strength H . For comparison, simulated curves with the unchanged resistivity parameters as obtained at 
fioH — 0.62 T are shown for hqH = 0.21 T (dashed lines). The transverse-resistivity parameters are found to be constant 
within the experimental magnetic field range. 




FIG. 8: Dependence of the longitudinal resistivity parameters 
po - P03 on the field strength noH for a (Ga,Mn)As(113)A 
layer with a Manganese concentration of 4.9%. The decrease 
of po with increasing /loH reflects the negative magnetoresis- 
tance. The lines are guides to the eye. 



the sample with 4.9% Mn are plotted as a function of the 
magnetic field. The longitudinal-resistivity parameters 
of the samples with 4.2% and 3.1% Mn showed a similar 
field dependence. In contrast, the transverse resistivities 
shown in Fig. [7] can be simulated with field- independent 
resistivity parameters pe — pn. In order to obtain a good 
fit of the experimental data, all parameters in Eq. (pij) 
with exception of pg and pn are required. The varia- 
tion of the lineshapes upon the field strength exclusively 
arises from the MA described by the anisotropy parame- 
ters from Eq. (I^IT) . Therefore, we mainly focused on the 
transverse resistivities when deriving MA parameters. 
Field-dependent resistivity parameters have also been re- 
ported by other groups for (113)A-oriented (Ga,Mn)A3^ 
and for (OOl)-oriented (Ga,Mn)A£24. In Ref.lHthe field 
dependence was studied up to 9 T. 

The microscopic origin of these findings is not clear yet. 
A (OOl)-oriented reference sample, grown at the same 
conditions as the (113)A-oriented layer with XMn — 4.9%, 
showed a similar field dependence of the longitudinal re- 
sistivities, indicating that the effect is not primarily re- 
lated to the substrate orientation. 



field-dependent longitudinal-resistivity parameters. In 
Fig. [HI the best-fit longitudinal-resistivity parameters of 
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IV. SUMMARY 

Starting from a continuum mechanical treatment of 
the lattice distortion in high-index epilayers, a general 
expression for the strain tensor e of (/i/iZ)-oriented layers 
was derived. The isotropic strain component Eh (and thus 
e) as well as the shear angle were related to the exper- 
imentally accessible quantities Ad/d and At. Applying 
the equations to the special case of (113)A orientation, eh 
and (/) could be experimentally determined for a series of 
(113)A-oriented (Ga,Mn)As layers using HRXRD. Based 
on symmetry considerations, analytical expressions for 
the free-energy density and the resistivity tensor were de- 
rived by means of series expansions in terms of the mag- 
netization components up to the fourth order, allowing 
for a phenomenological description of the MA and AMR, 
respectively. The anisotropy parameters were explicitly 
given as a function of the strain-tensor components. The 
expression for the resistivity tensor, deduced for mono- 
clinic and orthorhombic crystal symmetry, can be used to 
calculate the longitudinal and transverse resistivities for 
arbitrary current directions. In order to probe the MA 
and AMR of the (Ga,Mn)As samples by angle-dependent 
magnetotransport, expressions for the resistivities were 
derived for current direction along [332]. The measure- 
ments were performed at 4.2 K and revealed the presence 
of a strong uniaxial anisotropy 3^2771^^ along [001] which 
could be explained within our theoretical model by the 
explicit form of e. Further significant contributions to the 
MA were found to be B^im^, B^i.{m^+my), B^{m-n)'^ , 
and Bu{m • t)^. Whereas the transverse resistivity pa- 
rameters turned out to be nearly constant within the 
range of applied magnetic fields, the longitudinal resis- 
tivity parameters were found to strongly depend on the 
field strength. 



Appendix A: Partially Relaxed Layers 



instead of Eqs. ([5]) and 

From a reciprocal space map (RSM) around an asym- 
metric reflex (to 7^ n), At™ can be inferred, because it 
is the angle between the reciprocal lattice vector of layer 
GJ^j and substrate G^„, respectively, cf. Fig. [Q] Ej* can 
then be obtained from Eq. (|A3p . Assuming that the lat- 
tice parameter of the (real) substrate is known, the recip- 
rocal lattice of the substrate can serve as a reference; this 
allows an accurate measurement relative to the substrate 
without relying on the absolute angle scale of the diffrac- 
tomcter. The length of the layer's reciprocal lattice vec- 
tor I GJ„ I can be inferred from the RSM and consequently 
the lattice plane spacing of the layer = 27r/ |GJ„| is 



obtained. By inserting into Eq. (jA2l) . we find the 
value of d^. Because the virtual substrate is cubic, we 
obtain = d^Vh"^ + + P and with Eq. (jAip we find 
ai . Thus we can determine the degree of relaxation 



R 



ai 



(A4) 



If the relaxed layer is tilted with respect to the sub- 
strate, as it has been reported for relaxed (In,Ga)As lay- 
ers grown on (001)GaA£*, this tilt needs to be considered 
in the determination of ej^. The tilt angle ip = ^(Gj^, G^) 
can be inferred from a RSM around a symmetric reflec- 
tion {m — n), and the corrected angle At^" — Ar^ — i/j 
has to be inserted into Eq. ([6]) in order to obtain e^. 
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In this Appendix, we describe how the equations pre- 
sented in Sec. Ill AI can be used to characterize partially 
relaxed layers. Figure [3] illustrates that the layer in the 
partially relaxed state can be described as a layer which 
is commensurate with a virtual cubic substrate having a 
lattice constant Ov • The hydrostatic strain in the layer is 
then described by the parameter 



ai 



ai 



(AI) 



where ai is the relaxed lattice parameter of the layer. In 
order to describe partially relaxed layers, we have to re- 
place Eh by in all equations of Sec. Ill AI In particular, 
we obtain 



Ad 
~d 



Ar„ 



m 



■ cos T„ 



sm T„ 



(A2) 
(A3) 



FIG. 9: (color online) A partially relaxed layer can be thought 
of as a layer that pseudomorphically grows on a virtual cubic 
substrate with a lattice constant Ov different from the true 
substrate lattice constant as- Lattice planes m are parallel 
in real and virtual substrate. The relative inclination Ar™ 
of the plane m in the layer with respect to m in the virtual 
substrate is the same as the inclination with respect to m in 
the real substrate. 

If the layer was tilted with respect to the substrate the re- 
ciprocal lattice vectors and G\ of a symmetric reflection 
would not be parallel. Note that in this schematic the cubes 
do not necessarily represent the cubic unit cells. 



Appendix B: Resistivity tensor for monoclinic and 
orthorhombic symmetry 

We derived the resistivity tensors for monoclinic and 
orthorhombic symmetry up to the fourth order in m; 
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thereby we made use of von Neumann's principle as de- 
scribed in Ref. M- In Fig. [JUl we show that for (110)- 
oriented substrates the layer exhibits orthorhombic sym- 
metry. For this case the generating matrices are 



and 




[110]. [001] 



FIG. 10: (color online) Distortion of the fee eubie unit eell 
(blaek lines) when a layer with relaxed lattice constant ai is 
grown on a (110) plane of a substrate with lattice constant 
as < ai. The grey points denote the face centered atoms. 
For the zinkblende lattice the only symmetry elements are a 
twofold rotational axis around [001] and two planes; thus 
the point group is mm2 {C'2v)- The crystallographic unit cell 
of the body centered orthorhombic Bravais lattice is shown in 
red. Assuming pseudomorphic growth, the relations between 
the orthorhombic lattice parameters and the cubic lattice con- 
stant of the substrate are a = Os V2/2 and 6 — a,. 




(Bl) 




(B2) 



Note that here the matrix S5 has been adapted to the 
cubic frame of reference, where the symmetry operation 
is a reflection at the xy plane and not a reflection at the 
y plane as in the canonical representation for the matrix 
55. We obtain for the resistivity tensor 



Porthorhombic — ^tetragonal ~t~ 



orthorhombic 5 



(B3) 



where ptctragonai IS given by Eqs. (3), (4), and (5) in Ref. 
[Hand Aporthorhombic by Eq. (jBS]), respectively. 
For monoclinic symmetry, the only generating matrix is 
55 and we find 



clinic 



/^orthorhombic ~^ Apjuonochnici 



(B4) 



where Apnionociinic is given by Eq. (|B6p . The Greek let- 
ters in Eq. (jB5l) and (jB6l) are non-vanishing linear com- 
binations of the galvanomagnetic tensors, cf. Eq. (2) in 
Ref. [li 
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^Porthorhombic 



+ 



+ 



+ 



+ 



+ 



+ 








^imyniz liTUxniz ^zmxtriy 

52mxmymz 





€2mlml 





e-jm^mymz eTrnxmyUiz 



e-jrYixTriymz 
ermx-mlmz 




e^mlruz 

^2 



eQmxmy{ml + ml) 



(B5) 



monoclinic — 



+ 



+ 



+ 



+ 



+ 






-p2{mx + my) 
-Psrriz 



/32(mx +my) 


Parriz 



l^srriz 

-Pzmz 




75 TO J + 76 

75^2 + 76m2 








75TO^+76TO^ 75rn2+76m2 

jsm^mz + jgmymz Jwmzimx + my) 
liomzim^ + my) 'famymz + ^gm^mz 
^7mxmy 



-5i{ml+ml) 




-Saimlmy ■ 
—S^m^mz - 







5i{ml 



mxm 




drmlmz S^mlmz + Srmlm 




^rm^my 
^rm^my 
jnmzimx+my) 




S^m^mymz d^m^mymz 

Semlmz + drm^mz 
-Sem^mz - Srmlmz 



egm^ + egmy 







esm^ 



egm^ 



d^m^mymz 
-S^mxmymz 




€iiml 







eiim^mz + ti2mymz 



ei5mz{ml + my) 







eiQmlml eiom^m^ 

eiamlmy + e^mxm^ 



eiomlm'i 
eiomlml 




ei^mz{mx + my) e^mymz + ei2m^m2 e^mxmy + eum-^my 
eismlmy + e^m^m^ eism^m^ + eum%my ti(,mz{m% + m^) 



tnmxmymz + eigmxmymz 
eigmxmymzimx + my) 




eigmxmymz{mx + my) 
tnmxmymz + tmrnlmymz 




e2omxmymz{mx + my) 



(B6) 
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